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Abstract. A passive sub-millimeter precipitation retrievals algorithm is provided based on 
Microwave Humidity and Temperature Sounder (MWHTS) onboard the Chinese Feng Yun 3C 
(FY-3C) satellite. Using the validated global reference physical model NCEP/WRF/ 
VDISORT), NCEP data per 6 hours are downloaded to run the Weather Research and Forecast 
model WRF, and derive the typical precipitation data from the whole world. The precipitation 
retrieval algorithm can operate either on land or on seawater for global. To simply the 
calculation procedure and save the training time, principle component analysis (PCA) was 
adapted to filter out the redundancy caused by scanning angle and surface effects, as well as 
system noise. According to the comparison and validation combing with other precipitation 
sources, it is demonstrated that the retrievals are reliable for surface precipitation rate higher 
than 0.1 mm/h at 15km resolution. 
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1. Introduction 

Global monitoring of precipitation is important because of its significant human consequences. 
However, the multiplicity of hydrometeor types and their small-scale and large-scale spatial 
inhomogeneities make accurate measurements difficult. For example, rain gauge measurements are 
significantly impaired by wind, poor global coverage, and the non-uniformity of rain. Both ground- 
based radars and passive microwave satellite sensors sense precipitation aloft and are generally unable 
to discern how much of that precipitation evaporates before impact. Both are also sensitive to 
unknown local hydrometeor size, form, and vertical velocity distributions, as are simple single 
frequency radars on satellites. The resulting lack of adequate ground truth seriously complicates 
development and validation of global precipitation sensing methods. 

This paper addresses the following: 1) the consistency between histograms of millimeter-wave 
atmospheric radiances observed by satellite instruments FY-3C MWTHS [1-2] and those predicted by 
a mesoscale numerical weather prediction (NWP) model WRF 
(http://www2.mmm.ucar.edu/wrf/users/ ) in combination with a four-stream radiative transfer model 
(VDISORT). WRF was initialized by National Center for Environmental Prediction NCEP [3] 
analyses on a le grid approximately 6 h prior to each FY-3C MWHTS transit and employed the 
Goddard explicit cloud physics model. The scattering behavior of icy hydrometeors, including snow 
and graupel was considered. Retrievals of precipitation can have a coefficient of 0.61 now, and with 


the more and more observing data and validation data from radar and other techniques, this work is 
continuing and will have a better result. 


2. Instrument description 

MWHTS onboard FY-3C satellite is a fifteen-channel total power microwave radiometer in the range 
of 89GHz to 191GHz and has been launched in the September, 2013. Combining 13 horizontal 
polarized channels near 118.75GHz (which is the first operated to detect atmosphere on current 
operational polar orbiting satellite) and 183GHz, also window channels at 89GHz and 150GHz, 
precipitation content and path can be predicted. 

Data assimilation of MWHS is demonstrated to improve the analysis of numerical weather 
prediction (NWP) and play an important role in the ECMWF forecast mode [4]. Based on sounding 
channels operated between 89 to 191 GHz (horizontal and vertical polarization), the MWHTS 
operated in cross-track scanning mode and will be used for deriving atmospheric vertical distribution 
of temperature and humidity, rainfall, etc. The chosen observation channels are at 89 GHz (vertical 
polarization), 118.75 GHz (horizontal polarization, eight-channel), 150 GHz (vertical polarization), 
183.31 GHz (horizontal polarization, five-channel) [5]. 


Tablel Channel characteristics of MWHTS receivers 

No foo olarization E er n e — NEAT . se aoe 

“| “GHz” |? (MHz) | Pity |“ (KR) | RS | width 

l 89.0 V 1500 50 3—340 1.0 1.3 2.0° 
2 | 118.75+0.08 H 20 30 3—340 3.6 2.0 2.0° 
3 118.75+0.2 H 100 30 3—340 2.0 2.0 2.0° 
4 118.75+0.3 H 165 30 3—340 1.6 2.0 2.0° 
5 118.75+0.8 H 200 30 3—340 1.6 2.0 2.0° 
6 118.75+1.1 H 200 30 3—340 1.6 2.0 2.0° 
7 118.75+2.5 H 200 30 3—340 1.6 2.0 2.0° 
8 118.75+3.0 H 1000 30 3—340 1.0 2.0 2.0° 
9 118.75+5.0 H 2000 30 3—340 1.0 2.0 2.0° 
10 150.0 V 1500 50 3—340 1.0 1.3 1.1° 
11 183.3141 H 500 30 3—340 1.0 1.3 1.1° 
12 | 183.3141.8 H 700 30 3—340 1.0 1.3 1.1° 
13 183.3143 H 1000 30 3—340 1.0 1.3 1.19 
14 | 183.31+4.5 H 2000 30 3—340 1.0 1.3 1.1° 
15 183.3147 H 2000 30 3—340 1.0 1.3 1.1° 


“Note: In column 3, V/H means quasi V/H polarization 
3. Data collecting and processing 


3.1. Data collecting 
NCEP datasets(per 6 hours) from January, 1, 2014 to Oct, 31, 2014are chosen to run the Weather 
Research and Forecast model WRF, and derive the typical precipitation data from the whole world, 
especially for the extreme weather, such as typhoon, cyclone and rainfall. Bring out the incomplete 
cases and latitude values smaller than 0 degree and larger than 30 degree. The whole global datasets 
observed by FY-3C MWHTS from January, 1, 2014 to Oct, 31, 2014 are used to evaluated and 
analysis the validation of retrieval algorithm. 

Radiance data in Level-1B format from FY-3C MWHTS are employed for this study [6-7]. Both 
FY-3C has afternoon-configureured orbit. MWHTS is a self-calibrating total-power radiometer, which 
is a cross-track step-scan instrument and executes one complete revolution every 8/3 s period. The 
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MWHTS antenna systems have a nominal field of view (FOV) of 1.1 degree at the half-power points 
(covering a 16-km diameter footprint at nadir) and execute a cross-track scan with 98 Earth FOVs 
(within +53.35degreefrom nadir ), one cold space FOV, and one warm blackbody FOV per 8/3-s scan 
period. FOVs 1 and 98 (at +53.35degree from the nadir, respectively) are the outermost scan positions 
of the Earth views. The FOVs 49 and 50 (at +1.1° from nadir) straddle the nadir. The space view 
center is at 107.1°. The space views and internal warm target views are used to complete on-orbit two- 
point calibration for MWHTS. 

For FY-3C MWHTS, channel 9 and channel 15 are sensitive to precipitation. Combined with 
window channels distributed at 89GHz and 150GHz, people can retrieve precipitation rate using 
neural network estimators. 
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Figure. 1 Global brightness temperature values of 89 and 150GHz. 


3.2. Data processing 


3.2.1. Surface classification. MWHTS-observed footprints were classified as snow-free land, seawater, 
or sea ice and snow-covered land using the surface classification algorithm adapted from Chinnawat’s 
paper [8][9](Surussavadee, C 2006, 2008), where inputs include MWHTS channels land 10, land/sea 
flag, and surface temperature, which has been validated. In this paper, the authors focus on MWHTS 
footprints classified as snow-free land or seawater for global. 


3.2.2. Matching. Microwave humidity and temperature sounder (MWHTS) are onboard the Chinese 
FY-3Csatellite is a cross-track scanning radiometer with a swath width of 2650km. MWHTS has a 
nadir footprint size of 15 km, and 98 fields of view per scan line. The absolute visiting period is 5.5 
days. The limitation threshold of matching domain is 0.05degree for both altitude and longitude with 
20 mins. 

The atmospheric and surface conditions are used for producing global simulations of brightness 
temperatures that are measured by MWHTS. The vertical profiles of temperature, specific humidity, 
and pressure, the surface parameters of surface skin temperature, 2-m wind speed, and wind direction 
from the National Center for Environmental Prediction (NCEP) 6-h reanalysis data are used as input to 
WRF. The NCEP/WRF data field shave a horizontal resolution 0.05°X 0.05° and 29 vertical levels, 
with the highest vertical level located near 0.1hPa. 


3.2.3. Correction for small brightness temperature biases between NCEP/WRF/VDISORT and 
MWHTS. Considering the calibration accuracy and model error, bias correction of brightness 
temperature are in needed. According to the bias calculations and error corrections about cold and 
warm target, variable target, nonlinearity error and system random noise, the primary calibration 
results can be derived, which are shown in Figure. 2 and 3, which meet the requirements of design and 


development of MWHTS onboard FY-C satellite. To demonstrate that system requirements are 
appropriate and can be met, the paper provides the analysis of brightness temperatures through 
calibration process and profile retrievals using observing data of MWHTS. 


280 T T T T 240 


i \ | al 


220% 


270+ 


| + 200} 


1 1 1 1 1 1 1 i 
0 500 1000 1500 2000 2500 0 500 1000 1500 2000 2500 


Figure.2 Comparisons of brightness temperature between observation and simulation. 


Figure.3 Bias distribution of brightness temperature between observation and simulation. 


3.2.4. Omission of extreme values. The following values are omitted: 1) any brightness temperature 
for that footprint is less than 50 K or greater than 400 K, which is invalid, 2) the surface altitude is 
above 2 km for |latitude| < 600, or above 1.5 km for 600<|latitude| < 700, or above 0.5 km elsewhere, 
which could be covered by snow and are sensed more strongly and is called too-high, or 3) MWHTS 
channel 9 is less than 190 K, which implies that the atmosphere is so cold and potentially dry that 
precipitation is unlikely and that even the most opaque channel, i.e., 183+1 GHz, may sense the 
surface and yield false detections of precipitation, and is called too-cold. 


3.2.5. Principal component analysis (PCA). The method of principal component analysis is to extract 
only good signals that provide useful information about precipitation but are insensitive to most angle 
and surface effects, and other noises. Therefore, this method will filter out the redundancy caused by 


scanning angle and surface effects, as well as system noise. Also, the dimensions of data matrix are 
reduced, which is helpful to save time and make the retrievals easier. 


4. Retrievals and analysis 

ANN is essentially a nonlinear statistical regression between a set of predictors (in this case the 
observation vectors X) and a set of predictands (in this case profiles of atmospheric temperature Z). 
The ANN model consists 3 layers 1, 2, and 3, which represent the input layer, the hidden layer, and the 
output layer, respectively, as shown in Figure.4. 

The neurons of the input layer are represented by vector Xi (X1, X2X3, ++, XL), where L is the 
number of the input neurons. The neurons of the middle layer are represented by vector Yi (Y1, Y2, 
Y3, ***YM), where M is the number of the hidden neurons. The neurons of the output layer are 
represented by vector Zi (Z1, Z2, Z3), where Z is the rain detection and precipitation rate for different 
surface types and situations [10-12]. 
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Figure. 4 The schematic of artificial neural network. 


Estimates over land are usually slightly less accurate than those over sea, as expected, particularly 
at the lower rain rates where the surface remains visible. Using the artificial neural network retrieval 
algorithm of rain detection and precipitation presented above, and considering difference distributions 
of brightness temperatures as model inputs, as well as observing angles, the global map of rain 
detection can be retrieved as Figure.5 shows. 

In Table2, for most cases, the root mean square error (rms) errors are among the rain rate range, 
while for rain rate larger than 30mm/h, the rms error is less than the lower bound, which suggests good 
retrieval performance. While in typhoon region, the retrievals of rain rate are slightly worse and need 
more training data to improve the retrieving model. 
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(a) Global distributions of brightness temperatures on Jan20, 2015for (a) channel 1(b) channel 10 for 
FY-3C MWHTS 
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(c) Brightness temperatures of (d) rain detection in certain region 
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Figure.5 Brightness temperatures and rain detection results based on MWHTS on time 22:44Jan, 20, 
2015 to 00:30 Jan, 20, 2015. 


Table2 Rain detection of typhoon domain by FY-3C MWHTS between Jan,1 to Oct,31, 2014 


Location Ra ae = h) Land rms/K aoe 
0.1-1 0.88 0.84 

1-5 3.95 3.54 

Typhoon 5-10 7.11 6.78 
domain 10-30 22.1 17.21 
30-50 24.21 21.23 

50-65 34.2 32.4 


The correlation coefficients between rain rate detection and retrievals were found to be 
approximately 0.79 and 0.81 for land and sea, respectively. 


5. Conclusion 

For the recent work, the surface is classified the land mark as land and sea, and using the observing 
data to test and validate the accuracy of the rain rate. The work is not enough, and the authors are 
doing further improvement, like considering the surface covered by snow, ice and rainforest. Also, 
according to convection and stratiform, and classifying different rain rate from 0.1 to 50 as detail as 
possible to retrieve the rain rate. The radar data are also needed to validate the accuracy of rain 
detection. All of above will be described in future paper. 
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